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Abstract 

We study the robustness of the quintessence tracking scenario in the context of more general 
cosmological models that derive from high-energy physics. We consider the effects of inclusion 
of multiple scalar fields, corrections to the Hubble expansion law (such as those that arise in 
brane cosmological models), and potentials that decay with expansion of the Universe. We find 
that in a successful tracking quintessence model the average equation of state must remain nearly 
constant. Overall, the conditions for successful tracking become more complex in these more general 
settings. Tracking can become more fragile in presence of multiple scalar fields, and more stable 
when temperature dependent potentials are present. Interestingly though, most of the cases where 
tracking is disrupted are those in which the cosmological model is itself non-viable due to other 
constraints. In this sense tracking remains robust in models that are cosmologically viable. 
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I. INTRODUCTION 



Recent measurements of the cosmic microwave background (CMB) anisotropy power 
spectrum by the MAT Q, BOOMERANG @ and MAXIMA @ experiments suggest the 
Universe is nearly spatially flat, with a total density, fltoto, very close to unity. At the 
same time, there is ample evidence from observations, including CMB power spectrum HJ, 
galaxy clustering statistics ||], peculiar velocities [§] and the baryon mass fraction in clusters 
of galaxies || [7| that the density of the clumped (ie: baryonic and dark) matter in the 
Universe is substantially lower, being of the order of 30% of the critical value (fl ma ttero ~ 
0.3). Additionally the evidence from the spectral and photometric observations of Type la 
Supernovae |§ seem to suggest that the Universe is undergoing accelerated expansion at the 
present epoch (see e.g. || and references therein). 

One way of explaining this seemingly diverse set of observations is to postulate that a 
substantial proportion of the energy density of the Universe is in the form of a dark compo- 
nent that makes up the difference between the critical and matter energy densities, which is 
smooth on cosmological scales and which possesses a negative pressure. Various alternatives 
have been put forward as candidates for this dark component. One such candidate is a 
cosmological constant. This choice, however, involves an undesirable fine tuning problem, 
in that the ratio of the cosmological constant and the matter energy densities in the early 
universe need to be set to an infinitesimal value to ensure their near-coincidence at the 
present epoch. 

An alternative - and arguably more attractive - candidate is quintessence. Though little is 
known about the actual composition of quintessence, it has been shown that, should it exist, 
quintessence can in general be modeled as a scalar field rolling in a potential 0, an approach 
we adopt here. Additionally, it has been shown that some scalar field quintessence models 
have the appealing property of possessing attractor-type solutions for which the quintessence 
energy density closely tracks the energy density of the rest of the universe through most of its 
history [10|, |llj. The presence of such attractor-type solutions implies that their asymptotic 
behavior is largely independent of initial conditions. This allows the quintessence energy in 
the early universe to be comparable to that of the rest of the universe, thus providing the 
possibility of removing the fine-tuning problem that exists with the cosmological constant. 

Ultimately, any successful cosmological model must be based on a theory of high energy 
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physics, such as the string theory, M-theory or supergravity. These models are much more 
complex than is generally allowed for in usual studies of tracking. For tracking quintessence 
to be truly free from fine-tuning problems, the tracking phenomena must be robust in the 
more complicated and complete cosmological settings that derive from such high energy 
physics theories. Such scenarios typically include many scalar fields (some of which may not 
have potentials suitable for tracking), a Hubble law different from that of the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model at high energies and scalar field effective po- 
tentials that depend explicitly on the scale factor (or equivalently, temperature). Tracking 
would be of limited use as a solution to the fine-tuning problem if it was destroyed by the 
inclusion of such effects. 

The aim of this paper is to study the robustness of tracking in a more complete framework, 
by considering the effects of the three types of generalizations mentioned above. Though 
these additions make the resulting cosmological model more complex, it is important to 
consider their effects because in a realistic cosmology motivated by high-energy physics they 
are likely to be present, and because their inclusion leads to qualitatively new effects as well 
as new constraints. 

The structure of the paper is as follows. In section II we begin with the equations of 
motion for n scalar fields, with possible scale-factor dependent potentials as well as a non- 
FLRW expansion rate, and derive the corresponding fixed points. In section |TTT| we discuss 
how the tracking attractor-type solutions can be understood as a shadowing of instantaneous 



fixed points. Section [TV] contains our analysis of the stability of this shadowing and the 
independence of the attractors from the initial conditions. Finally section |V| contains our 
conclusions. 



II. EQUATIONS OF MOTION AND FIXED POINTS 

In a general cosmological scenario inspired by high-energy physics, one expects the pres- 
ence of a number of additional ingredients, among them: multiple scalar fields, generalized 
Hubble expansion laws and potentials that decay with expansion. Multiple scalar fields [4>u 
% — 1, . . . , n) arise naturally in theories of high energy physics. Generalized Hubble expan- 
sion laws (with H gen ) different from the standard FLRW Hubble law H FLRW can arise from 
corrections arising from the specific model one is considering. Two concrete examples being 
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the modified Hubble law appearing in a braneworld scenario [12|, 0, [15], and the effect of 
varying the strength of gravity. Explicit dependence of the potential on the scale factor can 
come about as an effective potential due to the interaction of <pi with another field, which 
has been 'integrated out', but with an energy density that decays with expansion - thus 
making the effective potential it induces for fa also decay with expansion. Here we shall 
employ potentials of the form Vi((f>i) = func(fa) exp (—kiN), where ki are constants and N 
is the logarithm of the scale factor. Note that ki > represents an explicit scale-factor, or 
equivalently, temperature, dependence in the potential V* of the field. This does not allow 
for direct coupling between tracking fields, but rather coupling between tracking fields and 
the field that has been implicitly integrated out. The choice of fcj = corresponds to the 
usual form of a scale factor-independent tracking potential. 

The evolution of such cosmological models is governed by the equations [||] 

fa + 3H gen fa + V ij(f>i = 
Vi{4>i) = func(fa) exp {-kiN) 



Hflrw = ^§ - (pBG + Ej Pj) 
«M^ = -3(l + WBG ) = -3 ( 7BG ) 

Where pbg is the energy density of the background which has a constant equation of state 
Ibg = 1 + wbg-, Pj is the energy density of the field 4>j and a dot denotes differentiation 
with respect to physical time. 

To study the possibility of tracking in such systems, we start by assuming H gen = H FLRW 
and in order to make contact with previous work [|14|| , we introduce the following change of 
variables 



y % JPtot ^/Ptot 



(2) 



The evolution equations then become 



t = + sJl^iV! + fa [Ej ((1 " k f) 2xj + | (x* + yf)) + lBG (l - Ej + V) 

% = - k iVi- sjlhxm + | w [Xj ((i - k i) ^) + k f Of + yf)) + ibg (i - Ej + v. 

= -3 (Ibg) 



(3) 



where A, = -I^i and k 2 = 8ttG 



N • 

The generalization to the case of H gen ^ Hflrw can be made in the following way. The 
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instantaneous effect of H gen on the evolution of the (xi, yi) variables can be included through 
a rescaling which involves replacing A, in Eq. (||D with 

Aj es A,%^. (4) 

However, this transformation does not encompass the effect of the modified expansion rate 
on the evolution of A*. Thus such an effect can qualitatively alter the nature of the attractor 



solution (see for example ||T5|). We shall return to this case in section [HI. 

In the following we shall also use the physically transparent set of variables and 7^, 
in terms of which the above evolution equations become 



1 dQ^j 



-37^ + 3 (-ibg (l - Ej %j) + Ej %jl<t>j) = -3 (% - 7) 



%■ = 3 {% - 1) ( 2 " \n ~ i] ; n = W 1 (5) 

where 7^ is defined as 

7*i = I 1 - "g I 7^ + J (6) 

and the equation of state for the field 0, is defined as 

70i = 7- ^ 

We also find it useful to define the weighted average of the equation of state 7, which is the 
rate of decay of the total energy density of the universe, in the form 

7 = E + ^bg (l - E Q 4>3 j • ( 8 ) 

We start by briefly discussing the special case of models with fixed Aj (i.e. r, = -L-Lgidi 
I ) . In that case, the above system of equations becomes an autonomous system with a 
set of true fixed points, which are obtained by solving the equations ^ = ^ = or 
™~ - = 0, Vi. We have calculated these fixed points for the systems ([D and (H) 



dN dN 

and the results are summarized in Tables | and 0. Included are the existence conditions 
as well as the required ranges of the A, and fcj. These are the generalizations of the fixed 



points given by (Tj|, to the case of models with n scalar fields and generalized temperature 



dependent potentials (with ^ 7^ 0). Note that in the single field case, in the hi = limit, 



points A through E correspond to the fixed points found in [LJ], while fixed points F and 
G are new. 
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For Vt^j ^ all the fixed points satisfy 7^ = 7, which implies (7 — 7bg) (l — J2j = 0. 
As can be seen from Tables | and [T| the fixed points that are most interesting for cosmological 
model building fall into two groups: A and B. These exist for fc; 7^ 6 and \ 7^ and can be 
represented by 



rr — /3 7 3_ 

A * V 2 A, 



y ( 2-7)(7-^) 



(9) 

^ = V 2" A,; 

or 

(1 - 1) - 1) = |n*A?, 7^ = 7 (10) 

which are only valid for 7 > y. Note that the effective equation of state is 7^ = 7, which 
is independent of the value of ki. Interestingly, the actual equation of state 7^ will decrease 
as ki increases so as to maintain 7^ = 7: 

For all points other than types A and 5, the equation of state of the field is fixed to be 
either 2 or |, which must be the same as the average equation of state 7 (and the equation 
of state of the background •jbg if ^bg > 0). Thus, from the point of view of the other fields, 
which see only the expansion rate, one can always absorb the fields at points other than A 
and B by a redefinition of the background Qbg — ^bg + J2je{fixed-y} ^<t>ji an d relabeling the 
B point by A if &bg > 0. No generality is lost because this redefinition of the background 
has no effect on the expansion rate, and it is only through the expansion rate that one field 
can affect another. In the following discussion we shall drop the tilde in Obg, an d by Qbg 
we shall mean the total background Q plus the contributions from fa at fixed points other 
than A and B. Thus in the light of above discussion, we shall only concentrate on the fixed 
points of type A and B and consider the following two possible scenarios: 

1. VLbg > 0: In this case each field fa is at the fixed point of type A. The average 
equation of state is the same as the background (7 = Jbg) an d the the fields can track 
the background energy density, making this point the most relevant for quintessence. 



The fl^i are in this case given by Eq. (10) to be 



, (7ag -g(l-ti 



Af 
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For this arrangement to exist as a fixed point one requires 

(l - k i)(lBG ~ f ) 

A 



< } =E3 V ° j \r 3J <i (13) 



Note that modestly larger values of the ki make point B less likely and point A more 
likely (because it is harder for the fields, for a given set of Aj, to come to dominate the 
total energy density) . Thus the presence of ki > generally makes tracking (point A) 
more robust. 

2. Cl BG = 0: In this case each field fa is at the fixed point of type B and we have 

n* = 3 ^" °)j T ~ (14) 



Now summing over i, one finds 



where a and \i are given by 



1 



A 



^-6(7=1 (15) 



i?"?^?"' * S ?lF^f- (16) 
Note that /i is analogous to A in the single-field case with ki = 0. One can make an 
analogy with parallel resistors in electrostatics, whereby the smallest A 2 dominates /z 2 . 



In the case of ki = 0, the results of previous work concerning assisted-inflation [16 
are recovered. Note also that because ki < 6, /i~ 2 > a > 0. However, the ensemble of 



fixed points does not exist for this entire range of values: they only exist if 7 G [0, 2] 
and G [0, 1]. One has 7 = fi 2 -^ 2 - < 2, which, in turn, fixes all of the fl^j. This 
arrangement does not exist as a fixed point for — 6a < 1 or ^ — 6a > 1 

It turns out that point A is more interesting for tracking quintessence models than B, 
as the latter is generally ruled out by a number of observational constraints, such as Big 
bang nucleosynthesis (BBN) and structure formation [0, [L7|]. However, point B has been 
studied in scenarios where the domination of the scalar fields is desirable, such as in the 



case of assisted inflation [1J|. Point A exists for 7 = •Jbg, while point B is only stable for 
7 < Ibg (recall we are considering fixed Aj (Tj = 1); the situation changes significantly when 
Aj varies, as will be shown in the next section). We also note that points A and B both 
require y < 7, Vz. Of course if y > 7 for some fa, that field becomes unstable but will 
generally have its corresponding FIm —> 0, thus making it harmless as a source of instability. 
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III. TRACKING BY SHADOWING AN INSTANTANEOUS FIXED POINT 



In this section we study the possibility of tracking in presence of variable Aj. Before doing 
this for the more general system (|), we begin by introducing the notion of tracking in the 
context of the simpler case of single scalar field models and briefly discuss how tracking 
might solve the fine-tuning problem, as was shown in |ll], |i~4]l . 



A. Tracking in models with constant Aj 

Consider the case with one scalar field, 0, and with the usual form of the potential given 
by k = and 7^ = 7^. Assuming A to be a constant, the fixed points become 'true' fixed 
points, given by 

(0 M , 7M ) = (^,Tbg) , A 2 > 3 lBG 



(^b,T*b) = (1,t)» A2<3 ^ 



(17) 



1BG- 
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Table II: Points available for Qbg = ; 
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7 = 2 Vj n^- (7^ - 2) = 



For large enough A, the system will be attracted to point A, and since asymptotically 
l<t> — Jbg, this offers a plausible solution to the fine-tuning problem. However, with a 
constant A, Vt^ is also constant, which implies that one can not have a significant contribution 
from quintessence to the present energy density (0.$ ~ 0.7) and at the same time satisfy 
the nucleosynthesis or structure formation bounds (7|, [17]] (fl^ ^ 0.15). To make the model 
compatible with observations, a variable A is required which has decreased from a large value 
in the early universe to order of a few today. Making A variable, however, means that the 
fixed points are no longer true fixed points, and the question becomes what determines the 
asymptotic dynamics of the system. An interesting (and useful) feature of the evolution 
Eqs. (|3|) and (|5D is that they only involve the value of A and not its derivatives. As a 
result, the rate of change of the phase space vector (^^,7^) at a given instant is the same 
as its corresponding value in the constant A case, with the same value of A. In this sense 
one may talk about fixed points at a given time, or instantaneous fixed points that would 
instantaneously act as attractors (or repelers) for the typical trajectories of the system. One 
could then imagine that there exist dynamical settings with slow enough changes in A such 
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that the trajectories shadow or track a moving instantaneous fixed point (A or B in this 
case). This is essentially the tracking scenario, which can in principle allow the quintessence 
energy density to follow the energy density of the rest of the universe through most of 
its history, and also eliminate dependence on initial conditions. The issue then becomes 
under what conditions does tracking behavior occur. It turns out that further conditions 
are required for the tracking solutions to exist and be stable. In the case of one scalar field 
with <C 1, it has been shown that there exist class of potentials for which tracking takes 



place provided the corresponding T satisfies T > 1 and is nearly a constant [|TI 



B. Tracking in more general settings 

In this section we extend the above analysis of tracking with a single scalar field and 
constant A to generalized settings with variable Aj's as well as 

1. multiple tracking fields 

2. scale-dependent potentials with ki ^ 

3. non-FLRW Hubble laws (expansion rate different from that of a FLRW universe) 

As our interest is primarily in tracking quintessence models, the emphasis in the following 
discussion will be on the shadowing of attractor point A in Table |. As in the case of models 
with a single scalar field, one would expect a stable fixed point to still act as an instantaneous 
type attractor. When the Aj's are allowed to change slowly, the usual trajectories shadow 
these points. This shadowing amounts to the field being attracted to a surface, defined here 
by rj = 1. From the form of the equation of motion for 7^, it is easy to see that = 1 is an 
'attractor surface' (for 7^ > y, which we shall assume), which yields a fixed value of 7^. It 
is interesting to note that this r« = 1 condition is precisely equivalent to Eq. (|10"|) , although 
the latter was derived as a fixed point for a constant A« case. We therefore focus on tracking 
that shadows point A and maintains a constant equation of state (though in principle other 
forms of tracking may be possible) and take 

r\ - -7T 1^,, ^ = 1 (18) 



3 (7* - I) (l 
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Figure 1: Schematic plot of the 'instantaneous fixed point' (solid dot) and the the shadowing point 
(open circle), showing their movement in the (7^,^) plane. For decreasing Aj (r$ > 1), the 
instantaneous fixed point evolves to larger Q^. The shadowing point shadows this point, riding on 
the r = 1 line. It is at a smaller 7^ such that it can stay on the r = 1 line as it increases. Eventually 
A will decrease to the point where <\> dominates (O^ ~ 1). 

which ensures 7^ = const and in turn implies 



% ^3 / _| W_|,_ (lg) 



For the system of Eqs. ([5]) to close we also demand 



dlnVLsi „d\xi\i rflnAi _ 
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For a FLRW expansion rate, one has 



dN A 

which gives 



3/ y»- k i) (20) 



din A; / h \ . 

-2^^ * -6 I % - j J (1 - r 4 ) ^ 3 (7 - %) 

which in turn yields an expression for 7^ in the form 

k- 7 — — 

7<Ai ~ — ~ const. (21) 

' 3 2r r l v ; 

Now for 7^ to be constant, as is required, either both 7 and Tj must be constants, or 
alternatively each needs to be arranged in such a way that the ratio appearing in ( |21|) is 
constant. But in general there is no a-priori reason for the latter and we shall therefore not 
consider this possibility further. Since 7 is the weighted average of all of the 7^, for tracking 
to occur for any field, we need either 7^ ~ const ~ 7 (which does admit the Q^j = 1 
case) or fl^j <C 1. For quintessence, the first case can be consistent with the BBN bound 



of < 0.15 [17J. However, if one is considering tracking fields in other scenarios, such as 
during inflation, ~ 1 can also be acceptable. 

We have given a schematic sketch of the tracking scenario in Figs. 1 and 2, corresponding 
to k — and k > respectively. In each case we have plotted both the instantaneous fixed 
point and the shadowing point depicted by a solid dot and an open circle respectively. Note 
that in each case the shadowing point falls on the r = 1 line, but at a shifted 7^, such that 
it can remain on the r = 1 line, as the system (A) evolves. Figure 2 illustrates the two 
significant effects of having k > 0: the region 7^ < k/3 is excluded and the distance between 
the 7^'s of the instantaneous fixed point and the shadowing point is narrower. 

In the case of models with a generalized expansion rate (with H gen 7^ Hflrw, as for 
example in models which include brane corrections or changes in the strength of gravity), 
the overall effect is to change the amount of friction the fields feel. More precisely, as was 
mentioned above, the form of the Eq. (|3|) implies that the attractor-type solutions exist as 
before, but their properties are determined by a new effective logarithmic slope \ defined by 
Eq. (H). The i — th field will then instantaneously behave as if it had this effective logarithmic 
slope. However, the evolution of Aj as a function of iV is not determined by rescaling Aj (N) 
by Eq. (^j, since the A» evolution equation does not transform in this way. As a result, when 
a correction to the expansion rate is present, the qualitative nature of the attractor can be 
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understood from the value (or range of applicable values) of \. The robustness of tracking 
then depends on the nature of the modified attractor of the field (pi, which is determined by 
the value of \. Of course by the time of nucleosynthesis, H gen ~ Hflrw, and the attractors 
will then be unmodified until the present time. 

Whether a given modification to the expansion rate helps or harms tracking depends on 
the details of cosmological scenario in question. For example, a larger expansion rate in the 
early universe may cause an attractor to be at a smaller value of the field today where the 
slope of the potential is larger, and the field less able to dominate the energy content of the 
universe - or it may turn out that at a smaller field value the slope of the potential is less 
and thus the field is more likely to become dominant. Additionally, the transition from the 
non-FLRW attractor to the FLRW attractor may be abrupt, if the field can not shadow the 
attractor during the transition, in which case tracking can be disrupted. Thus, the effect 
of modifications to the Hubble law on the robustness of tracking depends on the details of 
the cosmological model - but given the necessary details, the effect can be determined by 
analyzing the attractors resulting from Eq. fl4|) 



C. Tracking and independence from initial conditions 

In general, the attractor of field <pj is unique and independent of initial conditions. We 
will show this by constructing an equation involving and Peg as the only dynamical 
variables, which does not include the initial conditions of any of the fields (f>j. One can 
then imagine solving this equation to obtain 0j [/9bgL or equivalently, <pj [N\. In general one 
would expect this equation to possess a single, monotonically varying solution - although 
special cases where this is not the case can undoubtedly be constructed. The important 
point is that this solution is independent of the initial conditions of the fields. 

To see this, the first step is to find 7 as a function of the field values <pj. Using its 



definition and Eqs. (|i9|) and (|20|) for each field, one can find a quadratic equation for 7 as 



a function of A,-, r,-, kj with j = 1, n, in the form 



f 



E - 7 (l + E 4 fa + I) j + (^ BG + £ = (22) 
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Figure 2: Schematic plot showing the effect of making the field's potential scale factor dependent 
(V (<fi,N) oc exp (—kN)). The effect of k > is to set a minimum effective equation of state (that 
is, decay rate of p^) for the field and to make it larger for the attracting point. Thus the shadowing 
point shadows the 'instantaneous fixed point' more closely and comes to dominate more slowly. 
However, for k > 37 tracking is impossible and the field dies away (fl^ — ► 0). In this way k can be 
said to act as a 'throttle' for quintessence. 



As this is a quadratic equation for 7 G [0, 2], there may exist 0, 1, or 2 solutions for 7. The 
existence and stability of such solutions will be dealt with in the next section. The key point 
here is that solutions are manifestly independent of initial conditions. Combining Eqs. (|I"9D 



where 




(23) 



B 3 = (2T 



i) {ibg -!) + !■ 
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and yields 

„ f3\ (i-|)(7-») i 

i^i-^J 2r . _ l U# [<p k \ - — PBG . (24) 

One can also write an expression for 

= i — "1 — = ttts; = [<l>k] (25) 

2(2r i -l)(l-^J 

Thus one arrives at the following equation 

^ [0*] = 1 (26) 

This will generally be a very difficult expression to solve. However, it simplifies greatly when 
Q^i <C 1, because the fields are only coupled to each other through the effects each has on the 
Hubble expansion rate, and for Vt^ <C 1 they effectively decouple, reducing the complexity 
of the system. In principal, the solution (if it exists, is physical, and is stable - issues we 
shall treat below) yields <pi [N] , % = l,n, which are independent of initial conditions. Thus 
if the solution exists (i.e. is self-consistent) and is stable, then it is independent of the initial 
conditions. 

The above argument depends on all of the fields being 'tracking fields' - that is, Eqs. (|19D 
and fl2(iD being consistent with = 1 and thus 7^ = const. Of course, not all the fields in 
a multi-field model need necessarily possess a potential that is compatible with tracking, in 
which case those fields would then not track. In such a case the argument given above for 
the independence from initial conditions is no longer valid. It is then impossible to make a 
general prediction about the existence or initial condition-independence of attractors for the 
tracking fields. However, as the fields affect one another through altering the background 
expansion rate (or altering the value and rate of change of 7), one expects that if the VL^ 
of the non-tracking fields stay small, or if they do not cause 7 to vary rapidly, they will 
not affect the attractor solution of the tracking fields, and leave intact the above result 
concerning independence from initial conditions, for the tracking fields. 

IV. STABILITY OF TRACKING/SHADOWING SOLUTIONS 

We have shown above that the existence of instantaneous fixed points for the evolution 
equations (|5|) is mathematically consistent. However, for these points to give rise to tracking 
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behavior they must also be attractors, that is, the shadowing point, determined by Eqs. fll8|) 
and ( |2lD must be stable to perturbations. To find conditions for this, we shall obtain 
perturbation equations by perturbing the full evolution equations around the position of 
the shadowing point. From the nature of the eigenmodes of the perturbation equations the 
stability properties of the shadowing point can then be deduced. However, for models with 
n fields, the problem of determining the perturbation eigenmodes becomes one of finding 
the eigenvalues for a 2n x 2n matrix, or solving a polynomial of order 2n. Fortunately, an 
exact solution is not necessary in order to address several important issues. These include 
the stability in the nearly decoupled limit (f^ <C 1) and the determination of when the 
shadowing point will become unstable. Since the fields (pi only affect each other through 
the expansion rate, one would expect that for small <C 1 the system will behave as n 
decoupled systems, which turns out to be the case. The stability conditions for the decoupled 
fields correspond to the generalization of those found in [|14|], with the added explicit scale 
dependence of the potential or corrections to the Hubble law taken into account. For the 
coupled case the results are also fairly intuitive: from the the shadowing conditions, one 
can see that 7^ should be nearly constant, and equal to a nearly constant multiple of the 
overall average equation of state 7. Thus we have a key stability condition, namely that 7 
must be nearly constant. Now when VL BG dominates this condition is satisfied. However, if 
^4>i (r, — 1) become of order unity for any field </>j, even if each 7^ is constant, the proportion 
of the total energy density in the field increases rapidly, resulting in a rapid decease in 7. 
Similarly, as J2j ^<fj approaches unity, the condition 7 = ^bg — const (and thus 7^ ~ const) 
can no longer be maintained. As one would expect, this manifests itself as a growing mode 
of the perturbations, which amounts to an instability. This disrupts tracking as J2j ^<t>j ~^ 1- 
The important question is whether perturbations have growing modes if Ylj ^<f>j < 1? In the 
next section we find that the answer to this question is negative, so long as 7 and Tj remain 
nearly constant. 

A. The perturbation equations 

Assume for the moment that all fields are trackers - that is expressions (|T%| ) and ( fH| ) hold 
for each field individually. Now to study the perturbations, take 

n - 1 = ei, (2r, - 1) - 1) = -5 t + (7 - 1) (27) 
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with the perturbations €i,5i <C 1. The resulting evolution equations for these perturbations 
become 



dN 



dSj 
dN 



-6,3(7-1) ((2-|)- 



- k i 
2r,-i 



(28) 



I dry 
~ dN 



Sj) dr, 

2r,-l dN' 



where the perturbation in 7 due to the perturbation of each field is given by 

A7=|E,^ 

-ej2 ( 7 bg - k f] 



S 3 



7-1 



■■—if 



(29) 



This is the source of the coupling of the perturbations. Furthermore, the rate of change of 
7 enters into the perturbation equation, which is given by 



3 



( k \ d 
(Tj - 1) 6 % - v (% - Ibg) + ~7Tt (%j ~ Ibg) 



+ 



d'jBG 

dN ' 



(30) 



3 j v , w ™ ■ 

Note that we are not presently taking 7^^ to be a constant - the reason for which will become 
clear when we consider the effect on stability by non-tracking fields by absorbing them into 
a redefinition of the background. The condition for the closing of the unperturbed system 
of equations is that the terms in the stability equation above that are not proportional to 
6j or 5j must vanish. Thus the following terms: 

rr ,\ rfln (2iW) d-f 

j 

1 



(31) 



Jy J " dKfa ' dN 

must be negligible compared to e 3 - or Sj. For the i — th field 0, to have a stable r\ 
tracking attractor, one requires that ~^f— ^ ^ 1 for all % as well as Ej &<t>j (F? -1)^1 
and ^jjf- <C 1. That is, if the and Si equations without these terms are stable (have only 
decaying modes, with negative real parts of eigenvalues) then the and Si will decay until 
they are of the order of these 'residual' quantities. Alternatively, one can say that the r\ = 1 
shadowing point is shifted by a small amount of the order these terms. Thus we proceed by 
assuming the terms that are not proportional to ej or Sj are much smaller than 6j and Sj, 
and analyze the resulting stability equation with these terms removed. This results in the 
following equations for e, and Si 



% = S^ + tidi + ejbj + Sjd 



3^3 



dSj 
dN 



(32) 
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with 

n , — 3 (o I ^7 2r fc- 3 a k k r fc -i 

«fc — 2 V ^2r fc -l ^ 3 2T k -l 



b k = -3Q # [ lBG -!f- ^Jt 



2V k -l 



and 



2r fe -i 



(33) 



2r fc -i 



A 7 = E e A' + S 3 d 3 

j 

In general, finding the eigenmodes of the above perturbation equations is equivalent to 
diagonalizing a 2n x 2n matrix, or finding the roots of a polynomial of order 2n given by 

_ ^ (3b k + d k c k 

^/3(/3-a fe )-|c fc - ^ 



Finding explicit expressions for the eigenmodes in terms of the coefficients ( |33| ) is neither 
feasible nor useful. Even for the simple 2-field case, one is faced with extremely messy 
expressions for the roots of a quartic equation. One, however, does not need to solve for all 
the /3s explicitly. It is sufficient to impose the condition that the real parts of all relevant 
/3s be negative. A mode is only relevant if its growth leads to the failure of tracking of the 
whole system. One would expect that some fields may not track, and at the same time not 
interfere with the tracking of other fields if they effectively decouple, having f2<^ < 1. In 
that case, these are irrelevant, or harmless instability modes. As a result, Re [(3] > may 
not necessarily signal an instability in the entire tracking system, as the fields only interact 
by altering the value of 7, and with fl^ <C 1 the failure of 0, to track would not effect 
tracking by the other fields. 

A possible way of studying the stability of the tracking is to recast the perturbation equa- 
tions p^ ) into a system of equations analogous to n coupled, damped harmonic oscillators 
(5i) in the form 

1 d 2 5i (-Oj) d8i 3 ^ 



(- Ci )dN 2 (-c t )dN 2 

If one identifies N with time, then the oscillator energy is a good measure of the deviation 
of the tracking fields from their shadowing condition (background) values. Thus there are n 
oscillators <5«, each of which can be excited in 2n possible modes. Note that due to its form, 
equation ( |35D is not in general derivable from a conservative Lagrangian. In particular, note 
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that the coefficients of the terms in F proportional to Sj do not in general form a symmetric 
matrix (and those proportional to 6j an antisymmetric matrix), as would be necessary for 
F to be derivable from an interaction potential of the form 



d dV 



dV 



Thus the equations of motion do not conserve 'energy' - that is, the perturbation amplitudes 
may decay or grow. However, by examining the implications of equation (|35|) , one can 
determine what will happen - and thereby determine the stability of the system. 

To do this, we shall employ the ansatz Si = J2i Aue@ lN and then solve mode-by-mode to 
obtain 



pf - aS - -a 



E 



.4 



ji 



—d 



(36) 



3 V T 

Because the right-hand side is independent of i, we can immediately determine the relative 
amplitudes for mode I of oscillators i and k: 



A, 



~Ci) 



0? 



aS - -Ci 



A 



ki 



PI 



akPi - -Ck 



i ± k 



(37) 



We begin by examining the decoupled system (F — > 0). For the case of models with n 
scalar fields, there are 2n eigenvalues /? ; 0± , given by 

1 



01 



o± 



ai ± \/af + 6q 



(38) 



From Eq. ( |37D it is clear that for the Si oscillator only the modes with frequencies (3\ are 
present. The zeroth order stability condition is the requirement that Re /? 8 0± < 0. Using 
(j38[), this can be seen to be satisfied if and only if both eij < and < 0. Solving these 
inequalities yields the following requirements for stability of uncoupled tracking 

■MO 



l<7, 



2T,- 



1 > — t- > 0. 

2+7-2% 



(39) 



Of course, for the zeroth order stability conditions to be relevant, we need J2j ^<j>j (Xj — 1) <C 
1 and J2j^(/>j ^ 1- Note that for ki 



0, the stability condition becomes 2Ti — 1 > 



2+7 



or 7^ < ±1 y- L , which for 7 w ^ BG agrees with that previously obtained in QTTJ for the case 



of a single tracking field. This is not surprising, as in the decoupled limit of J2j &><f>j ^ 1 
each field can be treated individually. This lower limit on 2Ti — 1 varies between 1 and |, 
where Tj = 1 is allowed unless 7 = 2 and Tj < 1 corresponds to a diminishing tracking field 
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('Jn l < ^ or & t° t> e a "useful quintessence" field, that is to have > (I\ > 1), the 
condition (|39|) is automatically satisfied. For fcj > and Tj > 1 tracking is still stable, and 
the lower limit on 21^ — 1 becomes smaller, taking a value closer to |. One can understand 
this intuitively as follows: ki > is another channel for the decay of (potential) energy in the 
field. Thus as less energy goes into kinetic energy, A* changes less rapidly, and the system is 
more able to maintain the r» = 1 condition. 

We now consider the effect of the coupling between the modes given by F ^ in Eq. (|36| ) . 
The coupling causes the mixing of the modes, but diagonalizing the coupled system, as 
we have seen, is not practical. This is not necessary, however, as by examining the un- 
diagonalized system we can draw the important conclusion that the shadowing condition is 



stable until J2k^<j>k ~ > 1- From Eq. (]37|) it is clear that for each oscillator Si, the modes 
with frequencies (3f are dominant for F <C 1, but there are now also 2 (n — 1) other modes 
present, with the frequencies of the other oscillators. The amplitudes of these other modes 
are suppressed by a factor of order J2k ^4>k (unless (3^ ~ $ , in which case An ~ Ajj^). 



One can see from Eq. ( |36|) that the diagonal mode An can be treated as an oscillator, but 
with a shifted mass 2 . We identify the zeroth order mode's mass 2 with the | (— q) > term 
in Eq. fl36|) . The effect of the interactions is then to provide a mass 2 shift of 

^ = E pZaaZl* - W) ■ ( 4 °) 

The condition for stability then becomes the requirement that the shifted mass not cause 
the mode to become growing, which implies Re [/%] < 0. If Am 2 { was real, then this would 
be the condition that the shifted mass 2 must be positive 

mLu = I i-Ci) + > 0. (41) 

However, in general Am| is complex since (3j are complex. The stability condition then 
becomes 



Re 



m totii 



—Im 
2 



m tot u 



2 



> 0. (42) 



This consideration does not significantly alter the story - that there is a stability radius of 
order unity around m\ = | (— c,-) in the complex plane. Thus in this way one can see that 
the system is guaranteed to be stable for any parameter values that ensure the inequality 
|Am|| < m\ to be roughly satisfied. 
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Thus far we have been considering the diagonal elements of a mass matrix. We now 
consider the off-diagonal elements which correspond to the off-diagonal modes, whose am- 
plitudes are controlled by Eq. (p7[). Until Y^ktyk ~ 1, these off-diagonal mode amplitudes 
can not become of the same order of the diagonal mode amplitudes, and therefore Eq. (|42|) 
remains satisfied. We have seen that for J2j ^<f>j < 1 the diagonal mode amplitudes do not 
grow, and the off-diagonal mode amplitudes are suppressed relative to these by a factor of 
order J2k The implication is that the shadowing conditions (|l~8l ) and (]2ll) will be stable 
so long as J2j Q<t>j < 1 • 

The question remains as to what happens if some fields in the model are not trackers, and 



therefore do not satisfy the conditions ( (Eg ) and (|2T|)? These fields must then be excluded 
from the set of perturbations {Si, e^} and therefore do not contribute to A7. Because they 
do not obey the tracking conditions, little can be said about such fields in general. They 
do, however, contribute to ■Jbg- If their Q is large enough such that they can cause an 
appreciable variation of 7, they will cause a deviation from the shadowing conditions as 
noted above for the residual terms (|3lT). 



B. Summary of stability conditions 

We have shown that tracking with the = 1 condition is possible for nearly constant 7. 
The constancy of 7 is the key to closing the zeroth order tracking equations, and therefore 
it is not surprising that a constant 7 is the key to the stability of the shadowing points, and 
that many of the stability conditions can ultimately be traced back to it. We note that even 
though in principle other forms of tracking may be possible under other conditions, the type 
of tracking considered here is that which has been commonly considered in the literature. 

Note the near constancy of 7 in turn requires nearly constant Tj together with 
J2i^4>i (Ft -1) < 1 and J2% ^<t>i < 1- We have seen that the perturbations e« and Si will 
decay until they are of order dln ^'~ 1 ' > , J2j ^4>j (Pj ~ 1) or J^j resulting in the r ; = 1 shad- 
owing point to be shifted by small amount of the order of these terms. The system is more 
stable in the small fl^ limit, since as Y^j ^4>j (J*j — 1) or J2j tyj approach unity, this in gen- 
eral causes 7 to vary rapidly, thus making the = 1 condition to fail and the tracking to 
cease. The weakly coupled case is analogous to n individual tracking fields, and thus it is 



not surprising that the constraints in the limit of ki — > reproduce those given in ||1 1|| , for 
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the single field case. Furthermore, in this limit 7 ss 7bg, and the stability of the shadowing 



condition for the field (pi, given by (O) and fl2~l|) require 



2 



H) 



% < 7, 2r 4 - 1 > -A — U- > 0. ( 43 ) 

3 " 2+7-2^- 

If the stability conditions are not met for the field (pi, then clearly it will not track. However, 
the instability is irrelevant if Tj < 1. This is because for Tj < 1, the fl^i decreases with 
expansion (as 7^ > 7) and the field <pi will effectively become irrelevant as — > 0. Thus 
in this case the tracking of the other fields will not be effected by the failure of the field (pi to 
track. On the other hand, if > 1, or if T; = 1 and is of order 1, then the instability in 
the field (pi will ultimately ruin tracking for all fields, as it will cause 7 to vary rapidly, and 
Ti — 1 can no longer be maintained in general for any of the other fields. Finally, we have 
seen that the presence of non-tracking fields will harm tracking if and only if they cause 7 
to vary rapidly. 

We note that our stability analysis was limited to the violation of the shadowing condi- 
tions given by fll8l) and (|2~T1) . A concern that we have not addressed is the possibility that 
the homogeneous distribution of the scalar field may be unstable to formation of spatial 



inhomogeneities - the so called 'Q-balls' ||18|| . This is an important issue that we hope to 
return to in future work. 



V. CONCLUSIONS 

We have argued that for tracking quintessence to truly solve the fine-tuning problem, 
the tracking phenomena must be robust in the more complicated and complete cosmological 
scenarios that derive from high energy physics. Such scenarios typically include multiple 
scalar fields (some of which may not have potentials suitable for tracking), a Hubble law 
at high energies different from that of the FLRW model, and scalar field potentials that 
depend explicitly on the scale factor. Tracking would be of limited use as a solution to the 
fine-tuning problem if it were easily ruined by such effects. 

In models with n scalar fields, we have found that = 1 tracking requires a nearly 
constant 7. If any one of the tracking fields develops a large enough (Ti — 1), or if 
J2j &<t>j - > 1 then in general 7 will vary rapidly, causing the tracking to fail for all fields. 
Also, if any of the non-tracking fields become a significant portion of the total energy density 
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then they can potentially cause 7 to vary rapidly. In general, therefore, = 1 tracking 
requires Tj to be nearly constant for the field fa to track, as well as f2<^ (Tj — 1) <C 1 for all 
fields, and for J2j^<f>j not to approach 1. 

In models with a modified expansion rate (such as those including the brane corrections 
or changes in the strength of gravity), we have found the attractor-type solutions can still 
exist, but their instantaneous nature is determined by the effective logarithmic slope Aj. 
When such a correction is present, the qualitative nature of the attractor can be understood 
from the value of Aj. However, the answer to the question of whether a given modification 
to the expansion rate helps or harms the robustness of tracking depends on the details 
of cosmological scenario in question - but given the necessary details, the effect can be 
determined by analyzing the attractors resulting from Eq. (|). 

In models with temperature dependent potentials (fcj > 0), we have found that this de- 
pendence can make tracking more robust. Roughly speaking, ki > 'throttles' quintessence 
- slowing the rate of increase of fl^, thus making fields with larger less harmful to both 
tracking stability, and to cosmological scenarios. On the other hand, very large values of ki 
spoil tracking for field fa. If Tj > 1, a value of < ki < 37 keeps (7 — 7^) small, and slows 
the increase of fi^, making tracking of the entire system more stable. For ki > 3j the field 
fa will fail to track, though it becomes an irrelevant field and tracking of the other fields is 
not harmed. 

In conclusion, we have found that the conditions for robustness of tracking becomes 
more complex in the more general settings that derive from high energy physics. Tracking 
can become more fragile with respect to some such added complexities to the model, as 
there are more constraints to be satisfied. This can for example be seen from the above 
analysis of stability in presence of multiple scalar fields. Other additions, such as temperature 
dependent potentials, on the other hand may make tracking more robust. Interestingly 
though, most of the cases where tracking is disrupted are those in which the cosmological 
model is itself non-viable due to other constraints. For example, although tracking is less 
robust for larger J2j tyj, this would generally make the model non- viable due to constraints 
such as nucleosynthesis and structure formation. Thus tracking seems to fair well in these 



23 



general settings, once we confine ourselves to viable cosmological models. 
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